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1. Introduction 

The aim of this work was to explore the applicability of the twisted-mass Wilson fermion 
formulation [|l|, ^ as described in the review by A. Shindler ^ for investigations of lattice QCD 
at non-zero temperature. The use of the staggered-fermion formulation has computational advan- 
tages [Q], but remains conceptually controversial On the other hand, the often used clover- 
improved Wilson fermion formulation requires to determine action and operator specific improve- 
ment coefficients. The twisted-mass formulation, combined with a tree-level Symanzik-improved 
gauge action ^, appears to be a challenging alternative for non-zero temperature lattice simu- 
lations, since it offers automatic ^ (a) -improvement by tuning the bare quark mass parameter only. 
It allows high-statistics simulations in the range of pion masses rriji > 270MeV. 

As a first step we had to characterize the phase structure of the model by locating the tran- 
sition/crossover lines and surfaces in the three-dimensional j3 — k: — /lo-space. The results of this 
study supporting a conjecture for the phase diagram by M. Creutz from the chiral perturbation the- 
ory point of view ^ ^ were already presented in Refs. [10, [ij]. Here we give an overview of 



the phase diagram but concentrating on the thermal transition surface. Moreover, we discuss a first 
feasibility study carried out at maximal twist for large quark mass. In the quenched case we are 
going to demonstrate that automatic iff (a) improvement also works in the finite-temperature case. 

2. The 3D phase diagram and the thermal transition 

The j8 — KT-phase diagram for two-flavour lattice QCD with clover-improved Wilson fermions 
has been thoroughly studied for small time-extent A^^ = 4,6 a few years ago by the CP-PACS 
collaboration ^2, 13]. A schematic view of the emerging phase structure is shown in the left panel 



of Fig. |l[ The cusp of the strong coupling Aoki phase (see 1 14, 15] and references therein) - the 
latter (in the infinite- volume limit) characterized by a non- vanishing expectation value (YtYs l'^ Y) 
indicating the spontaneous breakdown of a combined parity-flavour symmetry - seemed tightly 
connected with the thermal transition line K",()8). 

Here, we consider Wilson fermions with the additional twisted mass term /Iq W^Ys ¥■ For 
the gauge action the tree-level Symanzik-improved gauge action is employed. Having included 
the twisted mass term a more complicated 3D phase diagram has to be explored. For lattice sizes 



Nt: = S,Na = 16, we were able to show [ ]1 1[ ] that the Aoki phase ends somewhere inside the interval 
j8 = 3.0, . . . ,3.4 and, around jS = 3.4, becomes replaced by a region of metastabilities indicating 
a first order transition area (the shaded area in the right panel of Fig. [I]), a remnant of a transition 
known also in the zero-temperature case 1 16, 17, 1^, p^ ]. 

In what follows we are concentrating on the thermal transition seen at values jS > 3.65 and 
not too small /Iq (otherwise we are still running into the metastability region). Since the hopping 
parameter K and the twisted mass parameter /Iq are directly connected with the bare quark mass 

(2.1) 

we expect a cone-like structure of surfaces of equal physics around the critical chiral line K = 
JQ-(j8), Ho = 0. As a first step one can scan the phase diagram in a larger K"-range in order to see 
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Figure 1: Schematic view of the phase structure as seen in older investigations [ |13[ ] for a temporal lattice 
extent A^^ = 4, 6 (left) and as found in this work [ pi] ] with twisted-mass fermions in the P — K — /iQ— diagram 
forA^T = 8(right). 



how the thermal transition surface extends above K"c(jS). The result is shown in Figs. ^ For j3- 
values j3 = 3.4,3.45,3.65,3.75 from the steep rises of the Polyakov loop and from maxima of its 
susceptibility we observe very clear signals for a thermal transition in K. But additionally, for (j8 = 
3.75, /lo = 0.005), we see a tiny jc-interval around Kc = 0.166 where the Polyakov loop exhibits a 
comparably little maximum, which could have been easily overlooked. Thus, with rising K starting 
from values below jq. we pass through subsequent confinement-deconfinement, deconfinement- 
confinement transitions (or better crossovers) below and above , respectively, followed again by 
a confinement-deconfinement transition far above K^. The latter transition surface extends to the 
next fermion doubler region in the phase diagram. We have seen by additional j3 -scans at fixed 
k: > Kc{p) that the Creutz cone structure [^] that we are exploring is connected with the upper 
confinement-deconfinement transition by a phase transition surface bending upward in K at larger 
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Figure 2: JC-scans of the Polyakov loop (left) and Polyakov loop susceptibility (right) for various j3 -values 
(j3 = 3.4,3.45,3.65 for = 0.0068; jS = 3.75 for 0.005). Vertical lines mark K-ip ^ 3.75). 



Zooming into the region around K"c(j8) the Polyakov loop and its susceptibility behave as 
shown in Figs. ^ The maxima or shoulders of the Polyakov loop susceptibility shown in the 
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Figure 3: Zoom into the behaviour of the Polyakov loop (left) and its susceptibility (right) versus K for 
various j3 and /Xo = 0.005. Vertical lines from left to right mark the chiral critical values (^(jS) for j3 ~ 
3.9,3.8 and 3.75, respectively. 



right panel indicate smooth transitions or crossovers. For j3 = 3.75,/Xo = 0.005 this can be clearly 
seen in Fig. ^ where also the so-called pion norm has been considered. The Gaussian shape 
lines are fitted to highlight the position of the expected crossover. Note that at the given j3 and 
JQ-(jS) =0.166 the value /Iq = 0.005 can be related to a pion mass value m-n ~ 400MeV and a 
temperature T ~ 210MeV. This is close to values recently reported by the DIK collaboration [pO]. 
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Figure 4: Polyakov loop susceptibility (left) and pion norm (right) versus K both for j8 =3.75 and /Xq = 
0.005. 

We tried to figure out how far the crossover or transition cone extends in the ^Uo-direction. 
Although we collected O(IO^) HMC trajectories per point this turned out to be a quite difficult task, 
because of the weak and noisy signals seen in the plaquette and Polyakov loop susceptibilities, in 
the corresponding autocorrelation times as well as in the pion norm variable. From K"-scans for 
the Polyakov loop at j3 = 3.75 and various /iQ-values drawn in Fig. ^ we would like to conclude 
that the cone surface ends somewhere in the interval 0.014 < /Xq < 0.025. Corresponding fits of 
the ellipse shape distorted by lattice artifacts can be done with an expression for the quark or pion 
mass obtained at next-to-leading order in lattice chiral perturbation theory, but still have a quite 



large uncertainty []11[]. Therefore, we will not show them here. 

3. A first feasibility study at maximal twist 

So far we have not yet taken advantage of the expected ff{a) improvement. For fixed j3 and 
K = Kc{P) one would like to change /Xq in order to vary the physical quark or pion mass. Since the 
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Figure 5: Polyakov loop scans as function of K at various jUo for j3 — 3.75. 



statistical signals for the crossover turned out to be very noisy in this case, we instead decided to 
fix Ho and to vary jS and K = K"c(j8) accordingly. The values for K^(j8) can be estimated from the 
zero-temperature case (see e.g. [|6|, ^'). For Hq = 0.040 - which corresponds to a quite large pion 
mass value niji ~ IGeV - we have found the results shown in Figs. ^. The 'critical' value j3f = 3.88 
can be translated into Tc ~ 280MeV, which is again in the same ballpark in comparison with [^]. 
We conclude that this strategy to satisfy the requirements of an automatic ff^a) improvement seems 
to work. 
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Figure 6: Polyakov loop (left) and its susceptibility (right) versus j3 for maximal twist K = Kc{l5) at Hq 
0.040. 



4. Automatic i^{a) improvement at T 7^ 

Finally, we have checked that the automatic ff{a) improvement really holds in the finite- 
temperature case. In the quenched approximation we have computed the pseudoscalar screening 
mass for varying spatial and temporal linear lattice sizes while keeping the physical mass ratio of 
pseudoscalar and vector states and the physical temperature fixed. The results are plotted in Fig. 0. 
They demonstrate nicely a linear behaviour in the square of the lattice spacing a(j8). 

'We acknowledge the help of the ETM collaboration providing us also with data prior to publication. 
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(a/ro)= 



Figure 7: Pseudoscalar screening mass mpsro versus lattice spacing {a/ro)^ obtained for lattice sizes A^o- 
24, . . . ,32 and A^T = 6, . . . , 16 at fixed T/T^ = 0.655(5) and mps/my ~ 0.75. 



5. Conclusions 

We are convinced that with the present study the necessary prerequisites for a serious non-zero 
temperature analysis with twisted-mass Wilson fermions have been collected in a sufficient man- 
ner. The structure of the three-dimensional phase diagram has been explored in the physical range 
for the two-flavour case. Still it is difficult to locate the (pseudo-) critical behaviour or crossover 
at fixed (j8, k: = Kc{P)) along the direction of varying twisted-mass parameter Hq. Therefore, in 
a feasibility study, we have taken advantage of automatic &{a) improvement at fixed /Iq by pass- 
ing through the crossover phenomenon changing j8 and keeping close to the chiral critical hne 
(j8, (^.(jS)) for which we can rely on twisted-mass results at zero temperature. For the quenched 
case we have demonstrated that ff{a) improvement really works in the non-zero temperature set- 
ting . We are now in the position to start the determination of the critical temperature and of the 
equation of state with extrapolations to the limits of realistic light quark masses and to the con- 
tinuum. In order to reach smaller pion masses we continue our investigation with A^^ = 10, 12 on 
correspondingly larger spatial lattices. 
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